On the sensitivity of the dijet asymmetry to the physics of jet quenching 
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The appearance of monojets is among the most striking signature of jet quenching in the context of 
ultrarelativistic heavy-ion collisions. Experimentally, the disappearance of jets has been quantified 
by the ATLAS and CMS collaborations in terms of the dijet asymmetry observable Aj. While 
the experimental findings initially gave rise to claims that the measured Aj would challenge the 
radiative energy loss paradigm, the results of a systematic investigation of Aj in different models 
for the medium evolution and for the shower-medium interaction presented here suggest that the 
observed properties of Aj arise fairly generically and independent of specific model assumptions for 
a large class of reasonable models. This would imply that rather than posing a challenge to any 
particular model, the observable prompts the question what model dynamics is not compatible with 
the data. 



PACS numbers: 25.75.-q,25.75.Gz 



I. INTRODUCTION 



Hard partonic reactions in perturbative Quantum 
Chromodynamics (pQCD) lead (to leading order in the 
strong coupling ag) to highly virtual back-to-back par- 
ton configurations which subsequently evolve as parton 
showers, hadronize and are then experimentally observed 
as back-to-back jet events. One of the most striking man- 
ifestations of the interaction of hard partons with a QCD 
medium as created in ultrarelativistic heavy-ion colli- 
sions (URHIC) is the appearance of monojets in events 
in which one of the outgoing parton showers has lost such 
a high amount of energy into the medium that it is no 
longer identified by a jet finding algorithm. 

Such monojets have been observed in experiment in- 
directly via triggered two-particle correlations at RHIC 
energies e.g. by the STAR collaboration fl'|, but only re- 
cently directly on the level of reconstructed jets by the 
ATLAS and CMS coUaborations at LHC energies [2,[1]. 

The LHC experiments quantify the observation by fo- 
cusing on the imbalance of the energy reconstructed for 
two jets in opposite hemispheres. Specifically, in the case 
of the ATLAS experiment, calorimeter towers are clus- 
tered within a jet cone radius of i? < 0.4, a trigger con- 
dition on the energy Eti of the most energetic jet of 
Eti > 100 GeV is imposed and for valid trigger the en- 
ergy of the leading jet in the opposite hemisphere Et2 is 
measured with the additional condition Et2 > 25 GeV 
imposed. 

The asymmetry or dijet imbalance is then defined as 



Aj = 



Eti — Et2 
Eti + Et2 



(1) 



where Aj — indicates a perfect energy balance be- 
tween the two jets whereas Aj — 1 indicates a true mono- 
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jet event (note that the condition Et2 > 25 GeV prevents 
this situation in the ATLAS analysis). 

A number of models which were able to account for 
the measured imbalance were subsequently proposed. 
Such models range from parametric estimates [41 via 
schematic modelling of energy loss and jet finding I5| to 
sophisticated NLO pQCD modelling of jet evolution in a 
medium, however without detailed medium or jet-finding 
modelling Q and complete Monte Carlo (MC) modelling 
of a jet embedded into a bulk medium with realistic jet 
finding using the MARTINI code 0. The success of 
these models quickly established that the measured di- 
jet asymmetry does not pose a fundamental challenge to 
the radiative energy loss picture. 

However, at least an equally interesting question is 
what the information content of the Aj measurement is, 
i.e. what model properties are constrained by the observ- 
able. The aim of this work is to address this question by 
computing Aj under a number of different assumptions 
for both the parton-medium interaction and the evolu- 
tion of the bulk medium. 



II. METHODOLOGY 

A rigorous comparison of a jet quenching theory with 
e.g. the ATLAS measurement would in principle require 
the following steps: 1) generation of a bulk event using for 
instance hydrodynamics with event by event fluctuating 
initial conditions 2) generation of a hard event inside the 
same bulk event (since there is an observable correlation 
between medium fluctuation hotspots and hard vertices 
I^, the hard and soft physics can not be modelled inde- 
pendently in a rigorous treatment) and evolution of the 
perturbative shower inside the medium 3) computation 
of the back-reaction of the shower to the medium to ac- 
count for correlated soft physics inside the jet area (for 
instance elastic recoil of medium constituents) 4) simu- 
lating the detector response to turn particle information 
into calorimeter tower information 5) clustering and anal- 
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ysis of the calorimeter-level event as in the experimental 
procedure. 

The full set 1) - 5) corresponds to a numerically very 
demanding workflow, and consequently in the most am- 
bitious theory comparison so far [? ], only part of this 
has been included consistently using the bulk matter 
evolution code MUSIC without fluctuations and the in- 
medium shower evolution code MARTINI [7]. In the 
present work, we do not aim at a rigorous comparison 
with the data but rather at establishing the information 
content of an Aj measurement, therefore the following 
simplifications are assumed: 

• the event by event fluctuations of the medium are 
(for the purpose of embedding and propagating 
hard partons) replaced by an initial-state averaged 
hydrodynamical model 

• no back-reaction of the medium to energy- 
deposition by the shower is included, all energy 
transferred to the medium is assumed to be out- 
side the jet definition, i.e. 'lost' 

• a simple set of cuts on hadron species and Pt re- 
places the mapping of particles to calorimeter tow- 
ers 

• instead of the anti kx algorithm clustering the 
whole event, a simple cone drawn around the pri- 
mary parton direction with radius R < 0.4 is 
used to define the jet, the (by assumption) smooth 
medium background can be subtracted by defini- 
tion 

• medium (and calorimeter) fluctuations are then re- 
introduced as uncorrelated Gaussian fluctuations of 
the medium background energy inside the cone area 

While these assumptions have been made to mimick 
the behaviour of the published ATLAS baseline simu- 
lation the reader should be cautious to read any- 
thing more than qualitative trends into a comparison 
with data. Note however that qualitatively, the simpli- 
fied jet definition contains the same elements as the real 
calorimetric jet problem, e.g. hadron-species dependence 
of the calorimeter response or a momentum cut intro- 
duced by the deflection of low Pt charged hadrons by 
the magnetic field. 

Inside this framework, we test three different scenar- 
ios for the parton-medium interaction and three different 
hydrodynamical models for the medium evolution. Out 
of each model class, two of three model are known ruled 
out by some other set of data and are only used to test 
if Aj is able to discriminate the elements in which these 
models fail. 

Since the various test cases are easily distinguished by 
other observables (such as the single hadron suppression 
factor Raa)-, it is important to set up the calculation in 
such a way that only the information provided by Aj is 
probed. In particular, the dijet asymmetry is related to 



a conditional probability, i.e. given a trigger Eti > 100 
GeV, what is the probability to find an away side jet 
with Et2^ We're specifically interested in the informa- 
tion content of this conditional probability, not in the 
probability to find a trigger jet per event (which is re- 
lated to the suppression factor R''^^ of jets) in the first 
place. 

The dependence on the reduction of the trigger rate is 
removed as follows: The parameters of the only realistic 
model combination (YaJEM-DE in LHC hydrodynamics) 
are tuned such that the single inclusive charged hadron 
suppression factor Raa at 60 GeV is described (this im- 
plies typically partonic kinematics above 100 GeV as rel- 
evant for the dijet observation). The parameters of all 
other model combinations are then adjusted to create 
the same suppression of the jet > 100 GeV trigger rate, 
i.e. all scenarios by definition show the same suppression 
of the single jet rate. This procedure is very similar to 
what has been used to assess the information content of 
the dihadron correlation suppression factor Iaa in 0, [To| . 

III. THE MODEL 

The modelling thus consists of the MC simulation of a 
hard partonic back-to- back event, the choice of a medium 
evolution and embedding of the hard event into the 
medium, the propagation of the evolving parton show- 
ers through the medium and the clustering and binning 
according to the jet definition. 

A. The perturbative hard process 

In LO factorized pQCD, the production of two hard 
partons k, I is described by 

^^AB^kl+X ^^ij^kl 

, 2 J J — = V'a;i/iM(a;i,(3^)2^2/j7s(a;2,(3^) ;j — 

dp^dyidy2 ^ dt 

(2) 

where A and B stand for the colliding objects (protons 
or nuclei) and 2/1(2) is the rapidity of parton k{l). The 
distribution function of a parton type i in A at a mo- 
mentum fraction xi and a factorization scale Q ~ pt is 
/i/yi(a^ij Q^)- The distribution functions are different for 
free protons [ll|, llSl and nucleons in nuclei (l3l - [T5| . The 
fractional momenta of the colliding partons j, j are given 
by xi^2 — ^ (exp[±yi] + exp[±?;2])- Expressions for the 

pQCD subprocesses ^^^-^ — is,i,u) as a function of the 
parton Mandelstam variables s, t and u can be found e.g. 
in [16]. 

To account for various effects, including higher order 
pQCD radiation, transverse motion of partons in the nu- 
cleon (nuclear) wave function and effectively also the fact 
that hadronization is not a collinear process, the distri- 
bution is commonly folded with an intrinsic transverse 
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momentum kx with a Gaussian distribution, thus creat- 
ing a momentum imbalance between the two partons as 

PTi + PT2 = kT. 

These expressions are evaluated in a MC framework 
described in [9, lOj, resulting in an approximately back- 
to-back pair of hard partons with known parton types 
and momenta. 

B. Embedding into hydrodynamics 



factor contains a correction for relative motion between 
parton and fluid element, for instance it involves Lorentz 
contraction) . The only parameters which are not fixed by 
the hydrodynamical medium are the adjustible constants 
K (or Kd) which dial the overall strength of the medium 
modification. These are determined by a comparison to 
reference data or by the requirement of a constant jet 
rate suppression as described above. 



For the modelling of the medium created in URHIC 
at the LHC, we use a 2-|-ld ideal hydrodynamical 
framework which has successfully been tested against 
bulk spectra and already been used as a background to 
compute single hadron suppression at LHC energies . 
We contrast this with two hydrodynamical evolutions for 
RHIC energies, a 3+ld code \W\ and a 2+ld code [20^ 
which do clearly not reproduce bulk observables at LHC 
correctly. 

These frameworks differ substantially in the geometry 
of the simulated medium, for a comparison between 2+ld 
frameworks at RHIC and LHC see e.g. [l3|, for a compar- 
ison between 2-|-ld and 3-|-ld at RHIC see [2l| and are 
known to be discriminated by other high Pt observables, 
for instance the reaction plane angle dependence of sin- 
gle hadron suppression Raa (0) [211 by hard dihadron 
correlations [lO{ . 

For the actual embedding of parton showers into the 
medium, we assume that the distribution of vertices in 
the medium follows binary collision scaling as appropri- 
ate for a LO pQCD calculation. Thus, the probability 
density to find a vertex in the transverse plane is 



P{xo,yQ) = 



T^(ro + b/2)T^(ro-b/2) 

TAA{h) 



(3) 



where the thickness function is given in terms of Woods- 
Saxon distributions of the the nuclear density pAi^, z) as 
Ta{y) = J dzpA{i^, z) and T^^(b) is the standard nuclear 
overlap function TaaO^) = J c^^s T^(s)Tyi(s — b) for im- 
pact parameter b. We place the parton pair generated 
in the previous step at a probabilistically sampled vertex 
(2^0 J 2/0) sampled from this distribution with a random 
orientation with respect to the reaction plane. 

When embedding parton showers into the medium, hy- 
drodynamical parameters need to be related to transport 
coefficients which are probed by the shower. In the fol- 
lowing, this is always done by assuming the relevant co- 
efficients q (virtuality growth per unit pathlength) or D 
(mean energy loss per unit pathlength) are related to the 
energy density of the medium e via 



q[D]{0 = K[Kd] ■ 2 • e(C)3/^(coshp(C) - sinhp(C) cos^). 

. (4) 

Here C is the position coordinate of a parton in the 
medium, p is the transferse flow rapidity and V' the angle 
between parton propagation and flow direction (the last 



C. In-medium shower 

The in-medium shower for each parton given the path 
through the medium determined in the previous step is 
computed using the MC code YaJEM [2^^. YaJEM 
is a framework to compute medium-induced changes to 
the kinematics of the evolving parton shower where the 
medium is chacacterized in a phenomenological way by 
transport coefficients such as q oi D . The reliance on 
transport coefficients means that YaJEM does not resolve 
the medium microscopically. It is based on the PYSHOW 
algorithm [25] and has an option to simulate coherent 
medium-induced radiation as well as an incoherent drag 
force (as characteristic e.g. for elastic collisions with the 
medium constituents). In the following, we use YaJEM 
to set up three scenarios: 

The default scenario (YaJEM-DE) is characterized by 
medium-induced radiation with a ~ 10% component of 
elastic energy loss 26:]. The minimum scale Qp down 
to which the shower is evolved in the medium is set by 
formation time arguments as 24 ] 



(5) 



where E is the parton energy and L the in-medium 
length. This creates a length and energy dependence 
of the medium interaction which is indicated by a large 
set of high Pt observables [271 . Here, the small elastic 
component is required by the low momentum strength 
of dihadron correlations which is overestimated without 
this contribution |26|. 

In an alternative scenario (YaJEM), we set the scale 
down to which the shower is evolved to a constant 
Qf = 0.7 GeV. This changes pathlength and energy de- 
pendence and leads to disagreement with the high Pt 
observables both at RHIC ^ and LHC [H]. 

For a second alternative scenario (YaJEM-E), we as- 
sume that the medium effect is only manifest via inco- 
herent drag and has no induced radiation component, i.e. 
q = 0. This changes observables like subjet fraction or jet 
shapes to a behaviour opposite to an induced radiation 
281 ] — while induced radiation widens a jet, a 



scenario 



pure drag force scenario narrows the jet shape (cf. RAD 
vs. DRAG in [1^). The pathlength dependence of such a 
purely incoherent scenario is also incompatible with data 
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D. Jet clustering and binning 

The YaJEM output in terms of hadron jets is then clus- 
tered to compute P(Ejf.t\E, i, x, y, </>), i.e. the probabihty 
distribution for a given parton i with energy E produced 
at a position (x, y) in the transverse plane propagating 
at an angle to be registered as a jet of energy Ej^t after 
passage through the medium. 

Ejet is defined as the sum of the energy of hadrons ful- 
filling the following conditions: a) at an angle of i? = 0.4 
or less with the jet axis b) the hadron species is one 
of 7r+, 7r~, 7r°, if +, if 7 c) Pt > 1 GeV (similar 

criteria have been suggested in the context of jet mea- 
surements by the STAR collaboration [305). To account 
for uncorrelatcd fluctuations of the medium background, 
a AE is sampled from a Gaussian distribution centered 
around zero with width 20 GeV and is added to the jet 
energies in the presence of a medium. In jsTj , it has been 
pointed out that uncorrelatcd background fluctuations 
have a potentially strong influence on the observed dijet 
asymmetry and hence need to be accounted for. 

In this way, the near and away side jet energies Eti and 
Et2 can be computed from the partonic event and com- 
pared with the experimental trigger conditions. Given a 
triggered event, Aj is then be evaluated and the proce- 
dure is repeated with the generation of the next partonic 
event. 



IV. RESULTS 

Before looking at the outcome of the calculation, let 
us briefly review what we might reasonably expect to see 
by a comparison with dihadron correlations. Given the 
same suppression of the trigger rate, triggered dihadron 
correlations are sensitive to the pathlength dependence 
of energy loss [2^ and clearly able to rule out incoherent 
scenarios with linear pathlength dependence. They also 
have tomographic capability, i.e. they are able to distin- 
guish two different hydrodynamical background models 
(lOj . At the heart of these is a phenomenon commonly 
called surface bias — the distribution of vertices leading 
to a triggered event is very different from the unbiased 
binary collision distribution Eq. Q and does not peak 
in the medium center but towards the medium surface 
0. This in turn forces a longer than average path for 
the away side parton, and the precise balance between 
surface bias and pathlength dependence of energy loss is 
reflected in the strength of the dihadron correlation. 

In addition, the for low associate momentum bins, the 
dihadron correlation is able to trace the flow of energy to 
subleading particles and establish a lower bound on the 
rate with which energy is transferred into the bulk fluid 
medium [26]. In principle, all these phenomena can also 
be realized in dijet correlations. 



A. The baseline 

Even in the absence of a medium, there is some im- 
balance between the jet energy seen on the near and 
away side. The strength of the imbalance depends on 
the jet definition — energy-momentum conservation en- 
sures that the energy is balanced (up to intrinsic kx) 
when summing over the complete flow of energy and mo- 
mentum resulting from the original parton evolution, and 
hence any imbalance is caused by part of the flow of en- 
ergy and momentum going outside of the jet definition. 



2.76 ATeV 




FIG. 1: (Color online) The dijet asymmetry Ajet evaluated 
in the absence of any parton medium interaction assuming 
that all energy within a cone of radius R < 0.4 can be de- 
tected (ideal energy flow, green), applying the jet definition 
described in the text, but without bulk medium background 
fluctuations (pp), and in-medium baseline without final state 
interaction, but including uncorrelated Gaussian fiuctuations 
of the medium background energy inside the jet area. 

As evident from Fig. [1] this is not an academic discus- 
sion — the result with angular cut only and without any 
cuts on Pt or hadron species is quite different from the 
result for the full in-medium baseline where the whole set 
of cuts and fluctuations has been applied. 

Note again that while the jet definition has been cho- 
sen to resemble the ATLAS baseline simulation results, 
there is no reason to think that it really corresponds to 
the ATLAS detector response and clustering procedure, 
however the implication is that the probability distribu- 
tion for losing energy out of the jet definition is similar. 

B. Dependence on shower- medium interaction 
scenario 

The dependence of the distribution of Aj on the as- 
sumed scenario for shower-medium interaction is shown 
in Fig. [5] Within statistical and systematic uncertainties, 
all three scenarios are indistinguishable and the width 
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FIG. 2: (Color online) The dijet asymmetry Ajet computed in 
the 2-1- Id hydrodynamical framework for LHC [3 for three 
different scenarios of the interaction of the evolving shower 
with the medium (see text). Shown is also ATLAS data 
for reference, note that however the jet definitions are not 
identical. 

of the distibution qualitatively agrees with the ATLAS 
data. Given that both the pathlength dependence and 
the energy-momentum flow inside the modified jet is very 
different in all three scenarios, this is a somewhat remark- 
able outcome. The dijet asymmetry observable is neither 
strongly sensitive to the fact that YaJEM and YaJEM-E 
have a linear pathlength dependence which is incompat- 
ible with the angular dependence of single hadron sup- 
pression [23 | nor to the difference between widening and 
focusing of the jet shape [2^. Only for large asymmetries 
Aj > 0.6 there is a hint of different model dynamics, in- 
dicating that a purely elastic energy loss scenario into 
the medium is somewhat less efficient in creating large 
asymmetries. 

The suppression of the trigger jet rate as compared to 
the vacuum case as found in these scenarios corresponds 
to a suppression factor of R''^^ = 0.48±0.02 in the trigger 
momentum range (i.e. for jets of 100 to about 150 GeV), 
in line with ATLAS results for Rcp [32,]. 



C. Dependence on the hydrodynamical evolution 

In Fig. 13] the dependence on the assumed hydrody- 
namial modelling is shown. Within statistical uncertain- 
ties, there is no such dependence visible despite the obvi- 
ous differences in geometry, evolution time and freeze-out 
conditions when going from LHC to RHIC energies. Note 
however that by definition the differences in mean density 
between RHIC and LHC are not probed here, since in all 
cases the suppression of the trigger rate is required to be 
a constant. This implies that the quenching power scales 
(K and Ko) are adjusted to a much higher value for 



FIG. 3: (Color online) The dijet asymmetry Ajet computed 
with YaJEM-DE in three different hydrodynamical evolution 
scenarios for the background. Shown is also ATLAS data [2| 
for reference, note that however the jet definitions are not 
identical. 

the same value of the energy density in the RHIC hydro- 
dynamical runs. If the same relation between medium 
density and transport coefficients would be assumed in 
all cases, the results would be very different, but note 
that such an overall normalization factor is much easier 
probed in single hadron suppression. 

The obvious similarity of the outcome leaves only the 
conclusion that Aj is not sensitive to the medium geom- 
etry, i.e. the observable has no tomographic capability. 

D. Dependence on other factors 

The dependence of Aj on two other factors is illus- 
trated by Fig. m If no fluctuations of the uncorrelated 
background medium inside the jet area are applied, then 
the shape of the Aj distribution is found to be less broad. 
Similarly, if the medium transport coefficients are set to 
50% of the values used in the previous investigations, 
the resulting shape is quite different. However, note that 
both procedures alter the rate suppression of triggered 
events and are therefore not on par with the dependen- 
cies investigated in the previous sections. 

These results demonstrate that the dijet asymmetry is 
sensitive to the overall strength of shower-medium inter- 
action and that likewise that the proper treatment of the 
uncorrelated background fluctuations is important. 



E. Dependence on jet cone radius 

One of the experimental handles to control the amount 
of background fluctuations is the jet cone radius, and 
hence the systematics of the calculation with R is an im- 
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FIG. 4: (Color online) The dijet asymmetry Aj^t computed 
with YaJEM, assuming no fluctuations of the background or 
assuming only 50% of the medium quenching power. Shown 
is also ATLAS data 2] for reference, note that however the 
jet definitions are not identical. 



portant crosscheck. Since the radius parameter controls 
the jet area and hence the total amount of background 
probed, the effect of uncorrelated background fluctua- 
tions should scale linearly with the radius parameter, i.e. 
clustering jets with R = 0.2 should reduce background 
fluctuations by 50% of the strength seen at i? = 0.4. In 
the following, we use this scaling assumption for the un- 
correlated background. 

However, while reducing the cone radius reduces back- 
ground fluctuations, at the same time the probability for 
energy to fall outside the jet definition is increased. As 
a result, the fluctuations in clustered jet energy given a 
certain parton energy increase as compared to the case 
of i? = 0.4 (in principle even in the absence of a medium, 
but certainly in a more pronounced way in a medium). 

The net result of a variation in cone size depends on 
how well these two opposing effects cancel, since Aj is in 
the end a measure of all jet energy fluctuations regardless 
of their origin. 

The results using YaJEM-DE in the LHC 2+ld hy- 
drodynamical evolution are shown in Fig. [5j To first 
approximation the result indicates that the decrease of 
background fluctuations with R is compensated by the in- 
crease in fluctuations of jet energy inside the cone. How- 
ever, this cancellation isn't perfect and leads to small 
diffferences at low values of Aj. 

For the same reason, also the suppression of the trigger 
rate is almost unchanged between R = 0.4 and R = 0.2. 
For the smaller cone radius, the calculation yields R''^^^ — 
0.46 ±0.02 for jets in the trigger range, i.e. from 100 GeV 
to about 150 GeV, in line with the ATLAS results f32| . 

Note that this cancellation of background against jet 
energy fluctuations inside the cone can not occur in the 
absence of background fluctuations. The fact that ex- 



FIG. 5: (Color online) The dijet asymmetry Ajet computed 
with YaJEM-DE for two different choices of the jet cone ra- 
dius, R = 0.4 (the default selection) and R = 0.2. 

perimentally similar Rcp and Aj is observed for both 
R = 0.4 and R = 0.2 is a strong indication that both 
types of fluctuations are present, strengthening the mes- 
sage of jsij that uncorrelated background fluctuations 
need to be properly understood. 

V. DISCUSSION 

The apparent insensitivity of the dijet asymmetry ob- 
servable Aj to even gross changes in shower-medium in- 
teraction scenario or medium geometry is certainly at 
first glance more than somewhat surprising. However, 
upon closer inspection, several points can be identified 
which may help to explain the findings. 

A. The effect of clustering 

Jet definitions as used in jet finding algorithms are de- 
signed to make meaningful comparison between pQCD 
calculations valid at large momentum scales Q and ob- 
served hadron showers in terms of e.g. calorimetric mea- 
surements possible. For this purpose, jet definitions are 
created to systematically remove dependence to physics 
at low Q where pQCD starts breaking down, i.e. by con- 
struction jet definitions are relatively blind to enhanced 
soft gluon emission or the dynamics of hadronization. 

However, the medium modifications of a shower take 
place precisely at scales of the medium temperature, 
which in itself is not far from the QCD fundamental 
scale, Q ^ T ^ ^qcd- The implication is that vac- 
uum jet finding is designed to minimize the sensitivity to 
physics at the medium scale, and as a result jets are fairly 
robust against final state interaction with the medium 
when compared to e.g. leading high Pt hadrons. 
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B. Lack of surface bias 

The relative robustness against medium modification 
translates directly into a lack of surface bias. This is 
illustrated by Fig. [S] which shows the probability distri- 
bution of the location of a vertex in the transverse plane 
given a triggered event. 



YaJEM, LHC 2+ld hydro 




x[fm] 



thermal momentum scale explains the robustness men- 
tioned above. However, as pointed out in j^], this ar- 
gument no longer applies to soft gluons in the thermal 
momentum range, as they are easily scattered to differ- 
ent angles. This behaviour is fairly generic: whenever a 
model can produce soft gluons in the shower, these soft 
gluons are allowed to re-interact with the medium and 
there is energy and momentum exchange between shower 
and medium, the soft gluon angles becomes randomized 
no matter what dynamics is assumed in detail (in YaJEM 
for instance this effect has been observed the distortion 
of low Pt jet shapes long before Aj was measured [28j). 
In essence, once a gluon in the thermal momentum range 
is embedded into a thermal medium and interacts with 
the medium, it is difficult to find any argument why the 
gluon should not thermalize within a very short time, 
and so soft gluon production generically equals energy 
transport out of cone. 

We thus end up in a situation in which hard gluon 
emission is too rare to make an effect (and not governed 
by shower-medium physics in any case since the momen- 
tum scales are too different), the dynamics of semi-soft 
gluon emission is largely erased by the clustering process 
into jets and the dynamics of soft gluons is generic and 
independent of specific model assumptions. This explains 
why any model which is capable of suppressing the rate of 
triggered jets expects also a very similar dijet imbalance. 



FIG. 6: (Color online) The distribution of vertices leading to 
triggered events (here —x denotes the direction of the trigger 
jet). 

While the peak of such distributions is usually shifted 
off-center by ^ 4 fm in hadron-triggered correlations, no 
surface bias is visible in the case of dijet correlations. In 
other words, while the strength of medium modification 
and the in-medium pathlength are correlated in dijet as 
well as dihadron correlations, the strength of the correla- 
tion is much weaker due to the relative robustness of jets 
vs. single hadrons. The presence of strong uncorrelated 
background fluctuations 0(10 — 20) GeV in the jet area 
(with the typical scale of energy transported out of cone 
of the same order Q) further dilutes the position - en- 
ergy correlation and effectively removes any surface bias. 
As a result, the observable is insensitive to pathlength 
dependence or medium geometry. 



C. Generic fate of soft gluons 

While the suppression of single hadron production can 
be achieved by fairly collinear gluon emission, the sup- 
pression of a jet rate requires the emission of large angle 
gluons to shuffle energy out of the jet cone. Applying this 
argument to medium-induced gluon production above a 



D. The role of background and jet definition 

The above arguments suggest that the overall strength 
of the medium quenching power is the only jet-medium 
interaction related parameter constrained by the mea- 
sured distribution of Aj. However, judging from Figs. [T] 
and m the outcome is very sensitive to both the precise 
cuts used to relate hadrons to calorimeter jets and to the 
assumed amount of background fluctuations (the latter 
point has also been made in [3l|). This leads to the sur- 
prising conclusion that in an approach which models the 
whole event, the correct modelling of the initial state fluc- 
tuations of the hydrodynamical medium and their evo- 
lution is more important than the modelling of shower- 
medium interactions. 

The dependence on the way hadrons are related to 
calorimeter jets may then explain the somewhat differ- 
ent shapes of the Aj as observed by various groups 0-01 ■ 
This should again be taken as a sign of caution not to 
read too much into the (dis-)agreement of the results pre- 
sented in this work (and in other works) with the data. 
But the very observation that the observable is more sen- 
sitive to technical details of jet finding and treatment of 
fiuctuations than to the physics it is supposed to measure 
is clearly a reason to worry. 
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VI. CONCLUSIONS 

Following the reasoning presented in the previous sec- 
tion, it seems that observing the emission of soft gluons 
at the medium momentum scale is problematic in princi- 
ple, as such gluons would thermalize under very general 
conditions. A more meaningful question is then how the 
energy and momentum carried by such gluons is trans- 
ported by the bulk medium, i.e. by hydrodynamical ex- 
citations. 

The physics of medium modifications to a parton 
shower is tlic^ii contained in semi-soft gluons with pr ^ 
few T which are sufficiently hard not to be defiected to 
large angles by re-interactions with the medium. How- 
ever, the physics at siich momentum scales is best not 
observed via clustering into jets, as the clustering al- 
gorithms are designed to suppress sensitivity of observ- 
ables to physics at these scales. As a consequence, trig- 



gerd hadron correlations have much greater tomorgraphic 
power. 

It would thus appear that fully reconstructed jets do 
not make the best observables to do high Pt tomogra- 
phy of the medium or to determine the precise nature of 
shower-medium interaction, and that approches to focus 
on this physics specifically by suitable designed n-particle 
correlation are more likely to answer such questions. 
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